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ABSTRACT: 

The dynamical behaviors of a discrete-time SIR epidemic model are analyzed in this paper. Existence 
and stability of disease -free and endemic equilibria are investigated. Basic reproduction number for the 
system is computed. Numerical simulations show that the system has complex and rich dynamics and 
can exhibit complex patterns, depending on the parameters. Also the bifurcation diagrams are 
presented. 
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I. INTRODUCTION 

Infectious diseases have tremendous influence on human life. Every year millions of human beings 
suffer or die of various infectious diseases. Controlling infectious diseases has been an increasingly complex 
issue in recent years. Mathematical models have been extensively used in the study of spread of epidemics and 
diseases. The models aim to capture the major factors that are responsible for the progression of diseases, and 
can forecast how diseases spread. The results obtained from these models are useful in predicting how infectious 
diseases develop and spread. Kermack and McKendrick derived the simple classical SIR model in 1927. In the 
classical epidemiological SIR model, a population of total size N is divided into S susceptible numbers, / 

infective numbers, and R recovered numbers. 

II. FORMATTING YOUR PAPER 

Analytical and numerical studies of system of difference equations models in epidemiology have 
uncovered a number of interesting dynamical behaviors including stable steady states, oscillation, limit cycles 
and chaos. In recent years, several authors formulated and studied natural discrete approximation to the 
continuous epidemic models obtained by applying Euler's Method. In this paper, we consider the SIR epidemic 
model as a system of difference equations. In this model, the susceptible host population is assumed to follow 

the logistic growth with a specific growth rate constants. The force of infection is . The parameter a 

measures the inhibitory effect, y is the natural recovery rate of the infective individuals. Also t\ and A represent 
the per capita death rates of infective and recovered. Throughout this paper, we assume tj ^ 0. 

SCu + l) = r5(iO[l-5W]-T 

1 + oJ^nJ 

1 + ai [n) 
R(n -\rl) = 7/(n) + (1 - XjRln} 

We restrict our attention to the reduced model described by the system of following two equations. 

SC« + 13 = r5 W [1 - 5 W ] - ^^^^ (2] 

+ 1> = \^ ^ + (i - 0^ + r))/(Ti) 

1 + aS{n} 

A point x" is said to be a fixed point (or equilibrium point) of a map f if /(x*) = 3:*, It is important to develop 
qualitative or graphical methods to determine the behavior of the orbits near fixed points. We consider the 
equilibria of system. The system admits the Equilibrium points = (^^^o) (the disease free equilibrium 

, — —J where K = a -hy and A = — + au -h ray, (endemic 

equilibrium of the model). 
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III. DYNAMIC BEHAVIOR OF THE MODEL 

The purpose of this paper is to study the nonhnear dynamics of system. Mathematically, the stability of 
the equilibrium points starts with linearization of the system (1). The stability of these equilibrium points can be 
investigated by linearization. The Jacobian Matrix for the reduced system (2) is 



r — 2tS — 



ff5 



(3) 



Trace ;Cy./.ir3 ^r-2rS .^^l^^z^+JL. + _ ^nd 



Dct !(S,LR) ={r-2rS- 
Lemma 1 . Let 



(4) 



be the characteristic equation for a matrix defined by (3). Then the following statements are true: 

• If every root of equation (4) has absolute value less than one, then the equilibrium point of the system 
(1) is locally asymptotically stable and equilibrium point is called a sink. 

• If every root of equation (4) has absolute value greater than one, then the equilibrium point of the 
system (1) is a source. 

• The equilibrium point of system (1) is called hyperbolic if no root of equation (4) has absolute value 
equal to one. If there exists a root of equation (4) with absolute value equal to one, then the equilibrium 
point is called non-hyperbolic. 



Proposition 2. The disease free equilibrium point is a 



Smk if 1 < r < 3 and - E — - ] < S < — — - 

Source if r > 3 and : — r - y > pL> 2 ... - y 



Cr-l3 



Non hyperbolic if r = 1 and 



r-1 



Proof: 



-ir+ 2 
0 



(5) 



The Eigen values of the matrix /(Fp) are -^-i = — r -h 2, An = + (l — K). In view of Lemma, we see 



that, Fn, is a sink if 1 < r < 



3 and 



J ^ P 



r-L 



r-1 



£"0 is a source if r > 3 and 



— y>ii> Z + — y. Also flp is non-hyperbolic it r = 1 and '- . 

The basic reproduction number if p is defined as the average number of secondary infections that occur 
when one infective is introduced into a completely susceptible host population. is also called the basic 

reproduction ratio or basic reproductive rate. For our system = ^^^^^^ 0}"Qi+y3 reproduction number. 



r — 



0k{i +r-^)a +a) 



il-K-) 



(6) 



Trace = = — 



A^aK 
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2rK 



T — 



8K 



IV. NUMERICAL EXAMPLES 

In this section, we present some numerical simulations of system (1) to illustrate our results. The 
numerical examples substantiate the analytical results obtained in the previous section. Numerical study of 
nonlinear discrete dynamical systems gives an insight in to dynamical characteristics. The time plots for S(n); 
I(n); R(n), phase portraits and bifurcation diagrams for the system (1) are presented. Dynamic behaviors of the 
system (1) about the equilibrium points under different sets of parameter values are studied. 
Example 1 . Take r = 2.95../? = O.S.y = 0.15.. u = 0.00999. 

Since if p then the disease-free equlibrium is globally asymptotically stable for any period. 




Figure 1. Dynamical Behavior of Trivial Equilibrium points < 1 
Example 2. Take r = 2.7..ff = O.SS.r = 0.059.. u = 0.00999. 

Since if p > 1 then the endemic equilibrium is globally asymptotically stable for any period. 

20 1 1 1 1 1 1 
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r = 2.7; a=0.02; P = 0.35; y= 0.059; m = 0.00999 



Figure 2. Dynamical Behavior of Trivial Equilibrium points if p > 1 

Example 3. Take r = 2.95../? = 1.7.. 7 = 0.15, _^ = 0.00999. The system exhibits oscillatory behaviour and 
the phase portrait is reduced to a limit cycle. 
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t = 450; r = 2.95; a = 0.8; b=1.7; d = 0.15; e = 0.00999 
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Figure 3. Time plot of Interior Equilibrium Point 

Example 4. Take c- = 0.3..r = 2.95.^a = 0.009.. r = 0.05. A = 0.0099. Also 8 is given three different 
values S = 0.097, 0.063, 0.19 together with the other parameter such that ifp < 1. The infection dynamics is 
illustrated in the Figure-4. 
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V. BIFURCATION 

Bifurcation theory is the mathematical study of changes in the qualitative behavior of a given family 
and the solutions of a family of differential equations. Most commonly applied to the mathematical study of 
dynamical systems, a bifurcation occurs when a small smooth change made to the parameter values (the 
bifurcation parameters) of a system causes a sudden 'qualitative' or topological change in its behaviour. In 
general, the term bifurcation refers to the phenomenon of a system exhibiting new dynamical behavior as a 
parameter is varied. The numerical analysis of bifurcation problems is concerned with the stable, reliable and 
efficient computation of solutions to nonlinear problems. The technique of linearization used for stability 
analysis may fail at bifurcation points, for near such points the dynamical behaviors of the system may differ 
qualitatively from those of its linearized system. The Neimark Saker bifurcation for discrete dynamical systems 
corresponds to the Hopf bifurcation in the continuous case. We use the numerical simulations to display 
bifurcation diagram for susceptible population. 

1 
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Figure 5. Bifurication 

The following diagrams exhibit the change in the dynamical nature of Sin] for values of r chosen from different 

ranges in the bifurcation diagram. 
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Figure 6. Dynamics in different values of r 



In the present paper we examined the bifurcation of a mathematical model for the spread of an 
infectious disease. In order to eradicate the disease, the basic reproduction number ifp must be lowered than a 

threshold. The basic reproductive number RO determines the existence of the equilibrium. When /fp < 1, model 
(1) has a unique disease-free equilibrium. When if p > 1, model (1) has a disease-free equilibrium and a unique 
endemic equilibrium. 
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